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ABSTRACT 

We present the results of the first four (quasi-)simultaneous radio (ATCA), X-ray (Swift, RXTE), and y-ray (INTEGRAL) observations 
of the black hole candidate IGR J1709 1-3624, performed in February and March 201 1. The X-ray analysis shows that the source was 
in the hard state, and then it transited to a soft intermediate state. We study the correlated radio/X-ray behaviour of this source for the 
first time. The radio counterpart to IGR J17091-3624 was detected during all four observations with the ATCA. In the hard state, the 
radio spectrum is typical of optically thick synchrotron emission from a self-absorbed compact jet. In the soft intermediate state, the 
detection of optically thin synchrotron emission is probably due to a discrete ejection event associated with the state transition. The 
position of IGR J17091-3624 in the radio versus X-ray luminosity diagram (aka fundamental plane) is compatible with that of the 
other black hole sources for distances greater than 1 1 kpc. IGR J17091-3624 also appears as a new member of the few sources that 
show a strong quenching of radio emission after the state transition. Using the estimated luminosity at the spectral transition from the 
hard state, and for a typical mass of 10M Q , we estimate a distance to the source between ~11 and ~17 kpc, compatible with the radio 
behaviour of the source. 

Key words. Accretion, accretion discs; X-rays: binaries; Radio continuum: stars, Stars: individuals: IGR J17091-3624, GRS 
1915+105, GX 339-4, H1743-322 



1. Introduction 

Most of the ~40 Galactic black hole binaries (BHB) are transient 
sources that are detected in X-rays during bright, months-to-year 
long outbursts. During this period, their luminosity varies by 
several orders of magnitude. At the same time, the broad band 
(~0. 1-200 keV) X-ray spectra show pivoting from the canonical 
hard state (HS) to the canonical soft state (SS), through inter- 
mediate flavours of these last two. These states also have clear 
signatures in the temporal domain (0.1-1000 Hz), with the ap- 
pearance of different types (dubbed A, B, C) of quasi periodic os- 
cillations (QPO) in the HS, har d intermediate (HIMS), and soft 
intermediate (SIMS) states (e.g.lRemiilard & McChntockll2006t 
iHoman et al 1 120051: iRemillard et alj|2002l) . 

BHBs are also known to produce strong radio emission, 
either from a powerful compact jet in the HS or from rela- 
tivistic discrete ejections at t he transition from the HIMS to 
the SIMS le.g.. lFendenl2006h . In the HS, a strong correlation 
between the X-ray and radio luminosity has been found by 
ICorbel et alj d2003l) who point out a very strong coupling be- 
tween the compact jets and the inner accretion flow. For Cyg 
X-l, this jet has been detected above ~500 keV (iLaurent et alj 
1201 lb . The LR a dio=£x 6 ration between the radio and X-ray lu- 
minosities is thought to be universal (aka the "standard" pattern 
or GX 339-4 like behaviour), although it is based on only eight 
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Galactic sources (among which GX 339-4 and Cy g X- 1 show a 
large quenching of the ra dio emission in the SS; ICorbel et al.l 
120031: (Fender et al.ll2010h . Seven sources have recently been 
shown to define a new track in the radio versus X-ray lumi- 
nosity diagram by following th e relation £Radio=£y 4 (hereafter 
the H 1743-322 like behaviour; iFender et alJl2010EICoriat et al.l 
1201 ll) . The total number of BHBs simultaneously followed in a 
multi-wavelength way is, however, still low. It is, thus, not clear 
which of the two types of behaviour (GX 339-4 or H1743-322- 
like) is the most common, or what leads to these differences. 
Multi-wavelength observations of any new active target, over 
the widest possible range of luminosities are, therefore, impor- 
tant for better understanding the jet behaviour and assessing the 
links between the accretion and ejection processes. To this end, 
we conducted Australian Telescope Compact Array (ATCA) and 
Rossi X-ray Timing Explorer (RXTE) observations as soon as 
IGR J 1709 1-3624 was seen to enter a new outburst in 201 1. 

IGR J17091-3624 was discovered with the INTErnational 
Gamma-Ray Astrophysics Laboratory (INTEGRAL) on Apri l 
14, 2003 while entering outburst dKuulkers et alj 120031) . 
Observations with other facilities (e.g. RXTE, Swift) permit- 
ted it to be classified as a low-mass X-ray binary, which 
was probably hosting a black hole (e.g. iLutovinov et aT1l2005h 
ICapitanio et al] 20091). The renewal of activity in late January 
201 1 dKrimm et al] 120111 iKrimm & Kenneal 1201 1[) prompted 
several observing campaigns on the source. It was, in partic- 
ular, found that IGR J 1709 1-3624 has similarities with GRS 
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1915+105, since it showed types of X-ray variability reminis- 
cent of those of the so-called "p-heartbeat" a nd "/?" classes seen 
in GRS 1915+105 dAltamirano et alj|201 lallbl) . The specifics of 
the IGR J17091-3624's X-ray behaviour further justify a (si- 
multaneous) multi-wavelength study, to see, in particular, how it 
compares to the other ~17 BHBs with known radio-X-ray be- 
haviour. 

Early in the outburst, we reported the pres ence of radio emis- 
sion compatible with a compact jet in the HS (ICorbel et al.l201 ll) 
and a ~0.1 Hz QPO at the same time (IRodriguez et al. 1201 1 ). 
Here, we present an in-depth analysis of four (quasi) simultane- 
ous radio-X-y-ray observations. To complement the ATCA and 
RXTE data, we use INTEGRAL and Swift observations to en- 
hance the spectral coverage. 

Table 1. Radio fluxes of IGR J17091-3624 during the four 
ATCA observations. Errors are at the lcr level. 



2011 February 



Obs. 


^5.5GHz 


^9GHz 


a 




(mJy) 


(mJy) 




Al 


1.40±0.05 


1.24±0.06 


-0.25+0.12 


A2 


1.53+0.10 


1.57±0.10 


+0.05±0.19 


A3 


2.41±0.10 


1.13+0.10 


-1.54±0.20 


A4 


0.17±0.05 


<0.08 





2. Observations and data reduction 

Figure [TJ shows the 15-50 keV daily averaged Swift/B AT light 
curve of IGR J 1709 1-3624, with the various observations dis- 
cussed in this paper. The journal of these observations can be 
found online in Table [4] The observations are, hereafter, or- 
dered and labelled according to the instrument used (A stands 
for ATCA, I for INTEGRAL, S for Swift, and R for RXTE). 
The ATCA observations were made in two frequency bands (5.5 
and 9 GHz) simultane ously, with the upg raded Compact Array 
Broadband Backend dWilson et alj 1201 lb . The amplitude and 
band-pass calibrator was PKS 1934-638, and the antenna's gain 
and phase calibration, as well as the polarisation leakage, were 
derived from regular observations of the nearby ( ~ 3.2°aw ay) 
calibrator PMN 1714-336. See, e.g.. ICorbel etatl ( 12001 for 
more details on the standard data reduction. 

The RXTE and Swift data were reduced with HEASOFlQ 
v6 . 1(8, and the INTEGRAL data with 0SA v9 . <0 The X-ray 
spectral fitting software package XSPEC 12 . 6 . Qq was used to 
fit the energy and power density spectra (PDS) after converting 
the latter into XSPEC readable files. 

The Swift/XRT level 2 cleaned event files were obtained from 
window timing (WT) mode data with xrtpipeline keeping 
only the grade events. We then extracted the spectra from 
a 40x10 pixel region centred on the best source position. The 
background was estimated from a region of the same size at an 
off-axis position. The ancillary response file (arf) was estimated 
with xrtmkarf, and the last version (v. 12) of the WT redistri- 
bution matrix file (rmf) was used in the spectral fits. 

The ZATOG/ML/ IBIS/ISGRI source spec tra were obtained in 
a standard way (e.g. IRodriguez et a l. 2008a). We used the latest 
version of the rmf and arf files available from the INTEGRAL 
instrument calibration tree provided with OSA 9.0. In two ob- 
servations, the source is also within the field of view (fov) of 
the JEM-X monitors, and we extracted 16 channels spectra from 
JEM-X unit 2. 
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Fig.l. 15-50 keV Swift/BAT lightcurve of IGR J17091-3624. 
The vertical arrows represent the (mid) times of the pointed ob- 
servations. The width of the arrow's top is the length of the ob- 
servation (visible only for INTEGRAL). 

The Z source GX 349+2 is at -41' from IGR J17091-3624. 
The RXTE/PCA data are thus contaminated by the emission from 
this bright source at a level that is very difficult to estimate 
given the lack of strictly simultaneous data. The RXTE/PCA 
data were therefore o nly considered for timing analysis (e.g. 
IRodriguez et al.l 1201 ll) . Two to 60 keV PDS were extracted 
from Good Xenon mode date with POWSPEC vl.8 over the 
0.0078125-500 Hz range. The fits of these PDSs were restricted 
to the frequency range 0.0078125-30 Hz. 

3. Results 

3.1. Radio observations 

The analysis of the ATCA data shows there is a single radio 
source within the X-ray error circle, at a location of {a, 5) = 17 h 
09 m 07 s .61, -36° 24' 25".7 (J2000, lcr positional uncertainty of 
0.1"). This position is also consis tent with the posi tion of the 
suggested optical/NIR counterpart dTorres et al.l201 ll) . The mea- 
sured radio fluxes and spectral indices a (defined by F^ocv") are 
reported in Table Q] 

3.2. X and y-ray spectral analysis 

The Swift and INTEGRAL spectra Sl+Il, S2+I2, S3+I3, and 
S4+I4 were fitted simultaneously, and we retained the energy 
channels 0.7-7.5 keV for XRT, 20-300 keV for ISGRI, and, in 
Sl+Il and II', the 4-20 keV channels for JEM-X2. The best- 
fit models consist of an absorbed power law with a high en- 
ergy cut-off (tbabs*cutoffpl in XSPEC, hereafter CPL) for 
spectra Sl+Il, IT, and S2+I2. A disc component is required 
in S3 +13 and S4+I4, while the cut-off is not needed any more 
(tbabs*(diskbb+powerlaw), hereafter DPL). A normalisa- 
tion constant was also included and frozen to 1 for the XRT 
spectra and left free to vary for the other instruments, except 
for the IT fits, where it is frozen to 1 for ISGRI. The best-fit 
parameters we obtained are reported in Table [2 and the 'v-F v ' 
spectra of Sl+Il and S3 +13 are represented in Fig. [2] 

In XRBs, a cut-off power-law spectrum is usually inter- 
preted as the signature of inverse Comptonisation of soft seed 
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Table 2. Results of the spectral fits to the joint SwiftfXRT and INTEGRAL/ISGR1 (and JEM-X in II and II') spectra. Errors are at 
the 90% level. 



Spectra model Nh 


kT M /kT,-„j 


r/r E cl „/kT, 


2 

Xv 


Unabs. Flux* 


(xlO 22 cm- 2 ) 


(keV) 


(keV) 

1 An , n nn inn+lV 


(dot) 


0.5-10 keV 3-9 keV 20-200 keV 



II' 

S2+I2 

S3+I3 
S4+I4° 



ThC 
CPL 
ThC 
CPL 
ThC 
DPL 
DPL 



i,iH -0.05 

1.1 frozen 
1.1 frozen 

1.15 ±0.06 
1.13 ±0.05 
1.5 ±0.1 

1.16 ±0.05 



0.2 (frozen) 
0.2 (frozen) 

0.2 (frozen) 

1.2 ±0.2 
1 90+0.01 
x - -0 03 



0.8 ±0.3 
1.49 ±0.09 

0.8 ±0.2 
1.54 ±0.05 
2.1^ 

2.3 ±0.3 

2 r - 1 



i " -13 

85!« 

91+17 
74+22 

107!| 
27!™ 



1.00 (89) 
1.0 (34) 

1.0 (34) 

1.1 (137) 
1.2(137) 
1.3 (128) 

1.2 (409) 



8.5 
NA 
NA 
12.1 
11.7 
38.2 
31.6 



4.3 
6.5 
7.2 
6.0 
6.0 
11.4 
12.4 



17.2 
22.4 
22.3 
20.8 
25.9 
7.3 
2.1 



In units of 10 10 erg cm 2 s 1 . Fluxes are normalised to 5wi/f/XRT, except II' normalised to ISGRI 
Observations made more than a day apart. Norm, constant is high and poorly constrained (> 2). 




0.01 7-1 



1 10 100 

Energy (keV) 

Fig. 2. Spectra of obs. SI +11 (XRT/JEM-X 2/ISGRI, grey filled 
circles) and S3+I3 (XRT/ISGRI, dark triangles). The dashed 
lines represent the best-fit models. 

photons by a thermalised (i.e., with velocities following a 
Maxwellian distribution) popul ation of electrons. A thermal 
Comptonisation model (comptt iTitarchukl 1 1 9 9 4h provided a 
good fit to Sl+Il, II', and S2+I2 (ThC in Table©. We assumed 
a disc geometry in XSPEC and fixed the temperature of the seed 
photons. This model also permits avoiding the divergence of the 
power-law flux towards low energy since Comptonisation signif- 
icantly contributes above 3xkTj„j and decreases quickly below. 

3.3. X-ray timing analysis 

None of the RXTE observations shows the "p-heartbeat" type of 
variability. Since GX 349+2 contaminates significantly Obs. Rl 
and R2 (Sec. 2), we do not discuss the level of continuum vari- 
ability since we cannot attribute it to IGR J 1709 1-3624 with 
certainty. The IGR J17091-3624 PDSs, however, show QPOs 
that are not associat e d with GX 349+2 dRodriguez et al.ll201ll 
IShaposhnikovll2TjTlt IPahari et alj[20TTh . To obtain their 'true' 
amplitudes during Obs. Rl and R2, we have estimated the net 
count rate of IGR J 1709 1-3624 in the RXTE/PCA by using 
our spectral results (Table |2| and simulated a PCA spectrum in 
XSPEC. The amplitude of the QPO is then A net =A raw x 
where A law is the rms amplitude of the QPO obtained from the fit 
to the PDS, CR net the simulated count rate of IGR J17091-3624, 
and CRpca the observed PCA rate. No correction was applied to 



R3a,b,c,d since GX 349+2 was out of the PCA fov. The QPO 
parameters obtained from the analysis are reported in Table [3] 
and two representative PDSs are plotted in Fig. [3] These QPOs 
can be classified as type C QPOs, except maybe those found in 
observations R3 that could be transitory between C and B. 

Table 3. Parameters of the QPOs seen during Obs. Rl, R2, and 
R3. Errors are given at the 90% confidence level. 



PDS v Q po Q A 



(Hz) (=v/FWHM) (% rms) 



Rl 


°'° 81 -« 


4.3±2.2 


12 3 +ly 
-hi 


R2 


o.ioo!™ 


4.2+1.5 


R3a 


4.57 ± 0.08 


4.2 ± 1.1 


8.8 ± 1.2 


R3b 


4.57 ± 0.08 


4.5 ± 1.3 




R3c 


4.28 ± 0.06 


2.4 ± 0.3 


12 5 +05 


R3d 


4.66 ± 0.05 


5.8 ± 1.5 


8.0 ±0.6 



g 0.01 



— ''T+fH** 
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Fig. 3. PDSs of observations R2 (top panel) and R3d (bottom 
panel). 

4. Discussion 

In the first two observations, the spectral parameters (Table |2]i 
and presence of type C QPOs a r e typical of a BH in the HS (e.g. 
iRemillard & McClintockll2TM iHoman et al.ll2005h . During the 
third observation, the source clearly is in a softer state. The tran- 
sition is associated with a discrete radio ejection (see below) as 
ofte n seen in BHBs (the cle arest example being GRS 1915+105, 
e.g. lRodriguez et aT1l2008al) . This, and the presence of QPOs at 
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4-5 Hz, indicates that IGR J 1709 1-3624 is in a SIMS at that 
time. No temporal (RXTE) data are available for the fourth ob- 
servation. The spectra are indicative of either a SIMS or a pure 
SS. 

The value of Nh obtained from the spectral fits is slightly 
higher than the 2°x2° weight ed averaged absorpti on in this di- 
rection (N H =6.0xl0 21 cm" 2 ; iKalberla et alll2005l) . This is con- 
sistent with a source lying at a rather large distance in the Galaxy. 

The HS— >SI MS state transition occurred between MJDs 
55612 and 55616 dPahariet al.ll201ll) . Obs. 2 is the observation 
in the HS that is the closest to the peak of the hard X-ray outburst 
(Fig. [1} and to the state transition. It has a bolometric flux of 
Fboi.trans~4.lxl0" 9 erg cm" 2 s" 1 (ThC model). Most BHBs have 
£t ra ns.HS^ss£4% L Edd ,i_200keV, wh en assuming a power-law spec- 
trum with T=2 (lYu & Yan 20090. Although this value slightly 
underestimates the bolometric luminosity at the transition, we 
take it as a 'hard' lower limit on L t rans, which therefore allows 
us to estimate a lower limit on the distance to the source. Using 
Fboi.trans estimated above, we obtain a dist ance <in7f»i~ll kpc 
for a BH of ~1OM . With L t ran S =10%L E dd (lEsin et al.lll997h we 
obtain djnoqi~l7 kpc0. Using a different m ethod (based on the 
QPO frequency -photon index correlation), Paharietaj] d201 ll) 
also favour a large distance with a range of masses between 8- 
11.4M . 

The radio spectra of observations Al and A2 are consis- 
tent with being flat and they indicate of the pre s ence of self- 
absor bed compact jets (e.g. ICorbel e t al. 2004} iFender et alj 
12004 . Observations A3 and A4 are characteristic of optically 
thin synchrotron emission and can be attributed to the emis- 
sion of a discrete ejection associated with the hard to soft 
transition. The observation of yet another BHB showing dis- 
crete ejections at the hard to soft transition probably indicates 
a link between transitions and (discrete) ejections potentially 
common t o all microquasars. In XTE J 1550-564 and GRS 
1915+1 05. lRodriguez et al J J20031 l2008allbl) have suggested that 
the Comptonised component was the source of the ejected ma- 
terial. Here, the change in spectral shape and drop in the hard 
X-ray flux (Table |2| is compatible with this interpretation. 

In Fig. [4] we plot the levels of radio and X-ray emission for 
IGR J 1709 1-3624 in c omparison with a sample of represen- 
tative BHBs (data from ICorbel et alj 120031 l2008h ICoriat et alj 
1201 It iRushton et all l201o¥ " assuming the distances estimated 
above. The two radio points with the lowest X-ray luminosities 
correspond to the HS, while the other two correspond to the op- 
tically thin radio flare. At a distance of 1 1 kpc, the source be- 
haviour joins the track followed by H1743-322 (Fig. [4]). In this 
case (and at any lower distance), however, the break of the ra- 
dio vs X-ray relatio n occurs at too low l uminosities compared to 
the standard BHBs dFenderet al.ll2010h . At 17 kpc the position 
of IGR J 1709 1-3624 is consistent with the radio/X-ray correla- 
tion of most BHBs (within the dispersion), although we cannot 
assign it precisely to a specific track (i.e. the GX 339-4 or the 
H 1743-322 one, Fig.|4|. In any cases, a large distance (> 11 
kpc) to the source is favoured by the radio-X-ray behaviour. 

IGR J17 091-3624 has re c ently b ee n compared t o GR S 
1915+105 dAltamirano et alj 12011 alibi : iPahari et al.l [20lTh . 
However, even at the favoured distance of 17 kpc, our simul- 
taneous radio-X-ray observations show that IGR J 1709 1-3624 



4 With exception o f Cyg X-l (HM XB) and GRO J1655-40 (dis- 
cussed as an outlier bygu & Yan 2 009!) . 

3 Fit of non simultaneous Swift/INTEGRAL observations performed 
on MJDs 55611-55612 leads to F bol-trims ~3.5xlO" 9 erg cm" 2 s" 1 . This 
lower value would place IGR J 1709 1-3624 at an even farther distance. 



does not reach the bright levels of radio luminosity observed for 
GRS 1915+105 (Fig.©, and is more similar to standard BHBs, 
such as GX 339-4 or H1743-322. This would tend to indicate 
that the radio properties discussed here are not related to the fast 
flaring X-ray behaviour common to both GRS 1915+105 and 
IGR Jl 709 1-3624. 
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Fig. 4. Radio vs X-ray luminosity for several BHBs. 
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Table 4. Journal of the multi-wavelength observations of 
IGR J17091-3624 presented in this paper. 



Facility 


label 


Obs. Id 


MJD start 


Good time 






(when available) 


(d) 




ATCA 


Al 




55601.64 


3.1 h 




A2 




55605.73 


1.9 h 




A3 




55623.57 


2.2 h 




A4 




55641.85 


2.4 h 


RXTE 


Rl 


96103-01-01-00 


55601.15 


2486 s 




R2 


96103-01-02-00 


55604.15 


5668 s 




R3a 


96420-01-01-010 


55623.05 


6755 s 




R3b 


96420-01-01-01 


55623.32 


4695 s 




R3c 


96420-01-01-020 


55623.49 


16263 s 




R3d 


96420-01-01-02 


55623.82 


7834 s 


Swift 


SI 


00031921005 


55601.05 


1457 s 




S2 


00031921009 


55605.13 


2167 s 




S3 


00031921019 


55623.45 


686 s 




S4 


00031921030 


55640.45 


2352 s 


INTEGRAL 


11 


Rev 1016* 


55599.49 


61080 s 




Il't 


Rev 1017 


55602.49 


131000 s 




12 


Rev 1018 


55605.48 


7818 s 




13 


Rev 1024 


55623.77 


6417 s 




14 


Rev 1030 


55641.98 


8583 s 



* INTEGRAL revolution number. 

1 II' is intermediate (in time) between Obs. II and Obs. 12. 



